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We study equilibrium properties of Bose-Condensed gases in a one- dimensional 
(ID) optical lattice at finite temperatures. We assume that an additional 
harmonic confinement is highly anisotropic, in which the confinement in the 
radial directions is much tighter than in the axial direction. We derive a 
quasi-lD model of the Gross- Pitaeavskii equation and the Bogoliubov equa- 
tions, and numerically solve these equations to obtain the condensate fraction 
as a function of the temperature. 

PACS numbers: 05.70 Ln, 05.70 Jk, 64. 



1. Introduction 

Recently ultracold atoms confined in optical lattices have attracted 
attention both theoretically^ - ^ and experimentally!^ - ^ In particular, the 
dipole-mode frequency^ and damping?' were measured by the direct obser- 
vation of an oscillating atomic current in a ID optical lattice at finite tem- 
peratures. Several papers have discussed finite-temperature properties of 
ultracold atoms in a ID optical lattice in addition to a highly anisotropic 
harmonic confinement, in which the radial confinement is much tighter than 
the axial confinement. Most theoretical studies have concentrated on the 
first Bloch band using the Bose-Hubbard model, and have ignored the effect 
of radial excitations P1EI This is only a good approximation when k^T <C -Br 
and k-gT <C h~uj±, where Br is the recoil energy that is roughly the high- 
est energy in the first Bloch band, and Huj± is the first excitation energy 
in the radial direction. In order to obtain more quantitative results that is 
applicable to the experiments of Refs. 6-8, however, it is important to con- 
sider thermal excitations in the radial direction, since the experiments do 
not always satisfy k-gT <C Br and k-gT <C h~LU±. 
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In this paper, we study equilibrium properties of Bose-condensed gases 
in a ID optical lattice at finite temperatures using Hartree-Fock-Bogoliubov- 
Popov (HFB-popov) approximation with taking into account the effect of 
radial excitations. We will show that the radial excitations are important in 
the thermodynamic properties at finite temperatures. 



2. Quasi ID modeling of a trapped Bose gas 

Our system is described by the following Hamiltonian : 



H = J dr\$(r) 



2m 



V 2 + Kxt(r) 



tjj(r) + — ip^ (r)^ (r)tp(r)tp(i 



(1) 



where g 



is the coupling constant determined by the s-wave scatter- 



ing length a. The external potential V ey ± is given by V eyi t(r) = Vtrap(r) + 



V op (z),"l4rap(r) 



[wj_(x 2 + y 2 ) + u 2 z 2 ] , V op (z) = sE^cos 2 (kz), where 



oj±,u! z are the frequencies of the harmonic trap, k = ^ is fixed by the 



wavelength of the laser field creating the ID lattice, 



2m 



is the recoil 



energy, and s represents the intensity of the laser beam. In order to take 
into account the quasi ID situation u>j_ > w 2 , we expand the field operator 
in terms of the radial wave function : 



(2) 



where (p a (x,y) is the eigenfunction of the radial part of the single-particle 
Hamiltonian, 



4>a{x,y) = e a (f) a (x,y), 



(3) 



(4) 



Here a = (n x , n y ) is the index of the single-particle state with the eigenvalue 
e (n x ,n y ) = fr^±( n x + n y + 1). Inserting (J2J) into (^Q) and using © and @, we 
obtain 



which satisfy the orthonormality condition 

dxdy 4>* a {x,y)4>j3(x : y) = 5 aj3 . 



h 2 d 2 m 2 2 TA , 



+ E 



9aa'/3l3' 



dz^lippiipa', (5) 
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where the renormalized coupling constant is defined by 

flW/3/9' = 9 dxdy4>* a 4>*p<t)f3i<t) a i- 



(6) 



Following the procedure described in Ref. 9, we separate out the field 
operator ip a into the condensate wavefunction and the noncondensed part; 
ipa = {ipa} + ipa = &a + ipa- Within the HFB-Popov approximation, we 
obtain the generalized Gross-Pitaevskii (GP) equation, 



h 2 d 2 m o ,, 



a' A3' 



(7) 



The Popov approximation neglects the anomalous correlation lip'tp)^ From 
the numerical solutions of the GP equation Q using the trap frequencies 
given later, we find that |$ Q | 2 /l$o| 2 < 1(T 6 (a / 0), and thus we will 
henceforth approximate <J> Q ~ $>5 a o, where we have denoted the lowest radial 
state as a = = (0, 0). 

Taking the usual Bogoliubov transformations for the noncondensate, 



3 

we obtain the coupled Bogoliubov equations, 

L a u ja + Y ( 2 9a' n o + gaa'pp'npp') u ja > - g%,n v ja 

a' - 

L a Vja + Y ( 2 9a' n + 9aa> W^w) Vja> ~ 9a' n 0Uja 
a' - 

where, 

m 



(8) 



EjUj a , 



(9) 



-EjUja, (10) 





ti 2 d 2 


2m dz' 


no(z) 


= i*wr, 




= SW00> 







2 + y w 2 z2 + V op (z) +e a -fi, 



(11) 
(12) 
(13) 
(14) 



Sums over the repeated indices (3, (3' are implied in Eqs. @, and pi)|). These 
equations define the quasi-particle excitation energies Ej and the quasi- 
particle amplitudes Uj and Vj. Using the solutions of Eqs. © and (fTUl) . 
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one can obtain the noncondensate density from 

(UjaUj/3 + V ja Vjf3) N{Ej) + V ja Vjfl 



ja/3 



(15) 



where N(Ej) = 1/ [exp(PEj) - 1]. 



3. Result 

Our calculation procedure is summarized as follows. Eq. (JJJ) is first 
solved self-consistently for ju and $ neglecting the interaction terms. Once 
<1> is known, Ej, Uj a and ujg are obtained from Eqs. @ and H10|) with ftp pi 
set to zero. This is inserted into Eq. and the process is repeated until 
convergence is reached. At each step, we define / dz& 2 = N — N, where 
N = J dzn is the total number of noncondensate atoms. 

Throughout this paper we use the following parametersPm( 87 Rb) = 1.44x 10 
kg, uj z /2it = 9.0 Hz, uj±/2it = 92 Hz, scattering length a = 5.82 x 10~ 9 m 
and the wavelength of the optical lattice A=795 nm. We fix the total number 
of atoms as iV = 4 x 10 5 . 
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Fig. 1. The condensate fraction N c /N as a function of the temperature. 
"s=0, g=0" represents the ideal Bose gas result without a lattice potential 
given by N c /N = 1 - (T/T°) 3 , where T c ° =141 nK. s=0, 0.2, and 0.4 are 
calculated using the HFB-popov equation. 



In Fig. ^ we plot the condensate fraction N c /N as a function of the 
temperature for various values of the lattice depth s. It can be seen that each 
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Fig. 2. The dipole mode frequency as a function of the lattice depth s with 
a fixed temperature T = 40nK. 



line falls to zero at approximately 140 nK, which is close to the semiclassical 
prediction of the BEC transition temperature of an ideal Bose gas in a 3D 
harmonic trap T c °= 0.947i(u;j_u; z ) 

1/3^1/3 /£. B= 141 nK _ In contrast, T c of 
an ideal Bose gas in a ID harmonic trap is T c ~ fiuj z ks ^ n ^ =13.4 //K. This 
shows that one must explicitly take into account the excitations in the radial 
direction in order to obtain the correct thermodynamic behavior at finite 
temperatures. 

The highest T c is obtained for an ideal gas without a lattice potential. 
The transition temperature deceases with increasing the lattice depth. This 
is due to a reduction of the excitation energy with increasing the lattice 
depth. The energy shift can be seen from the dipole-mode frequency shift. 
In Fig. |2]we plot the dipole-mode frequency as a function of the lattice depth 
s with a fixed temperature T = 40nK. One can see that the frequency shift 
from s=0 to 0.2 is bigger than the shift from s=0.2 to 0.4. This is consistent 
with the fact that T c shift from s=0 to 0.2 is bigger than that from s=0.2 
to 0.4. 

Fig. El shows the effect of the radial excitations in the noncondensate 
fraction. The solid line shows the noncondensate fraction, while the dash 
line represent the result obtained by neglecting the radial excitations, one 
can clearly see that it is necessary to take into account the excitations in the 
radial direction at almost all temperatures. 
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Fig. 3. The temperature dependence of the noncondensate fraction for s=0.2. 
The solid line is the result including the radial excitations, while the dashed 
line is the result neglecting the radial excitations. 

4. Conclusion 

We have presented a detailed calculation of the condensate fraction in 
a ID optical lattice at finite temperatures using the HFB-Popov approxima- 
tion. We explicitly included the radial excitations. Although the condensate 
is well described in terms of the lowest radial mode, it is necessary to take 
into account the thermal excitations in the radial direction at almost all 
temperatures. The shift of T c with increasing the lattice depth is due to 
the shift of the low-energy excitations. In a future study, we will extend 
the present work to calculate damping rate of collective modes and study its 
dependence on the lattice depth s. 
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